Quantum-limited linewidth of a chaotic laser cavity 
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A random-matrix theory is presented for the linewidth of a laser cavity in which the radiation is 
scattered chaotically. The linewidth is enhanced above the Schawlow-Townes value by the Peter- 
mann factor K, due to the non-orthogonality of the cavity modes. The factor K is expressed in 
terms of a non-Hermitian random matrix and its distribution is calculated exactly for the case that 
the cavity is coupled to the outside via a small opening. The average of K is found to depend 
non-analytically on the area of the opening, and to greatly exceed the most probable value. 
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It has been known since the conception of the laser 
that vacuum fluctuations of the electromagnetic field ul- 
timately limit the narrowing of the emission spectrum 
by laser action. This quantum-limited linewidth, or 
Schawlow-Townes linewidth, 



5u = \V 2 /I 



(1) 



is proportional to the square of the decay rate T of the 
lasing cavity mode Q and inversely proportional to the 
output power / (in units of photons/s). Many years later 
it was realised that the fundamental limit is larger 
than Eq. (|l]) by a factor K that characterises the non- 
orthogonality of the cavity modes. This excess noise 
factor, or Petermann factor, has generated an exten- 
sive literature (see the recent papers and references 
therein) , both because of its fundamental significance and 
because of its practical importance. 

Theories of the enhanced linewidth usually factorise 
K = KiK r into a longitudinal and transverse factor, 
assuming that the cavity mode is separable into longi- 
tudinal and transverse modes. Since a longitudinal or 
transverse mode is essentially one-dimensional, that is a 
major simplification. Separability breaks down if the cav- 
ity has an irregular shape or contains randomly placed 
scatterers. In the language of dynamical systems, one 
crosses over from integrable to chaotic dynamics [ [To|| . 
Chaotic laser cavities have attracted much interest re- 
cently jnj, but not in connection with the quantum- 
limited linewidth. 

In this paper we present a general theory for the Pe- 
termann factor in a system with chaotic dynamics, and 
apply it to the simplest case of a chaotic cavity radiat- 
ing through a small opening. Chaotic systems require a 
statistical treatment, so we compute the probability dis- 
tribution of K in an ensemble of cavities with small varia- 
tions in shape and size. We find that the average of K — 1 
depends non-analytically oc TlnT^" 1 on the transmission 
probability T through the opening, so that it is beyond 
the reach of simple perturbation theory. The most prob- 
able value of K ~ 1 is oc T , hence it is paramctrically 
smaller than the average. 



The spectral statistics of chaotic systems is described 
by random-matrix theory pO| , p^ |. We begin by reformu- 
lating the existing theories for the Petermann factor [^]|| 
in the framework of random-matrix theory. Modes of a 
closed cavity, in the absence of absorption or amplifica- 
tion, are eigenvalues of a Hermitian operator Hq. For a 
chaotic cavity, Hq can be modelled by an M x M Her- 
mitian matrix with independent Gaussian distributed el- 
ements. (The limit M — > oo at fixed spacing A of the 
modes is taken at the end of the calculation.) The ma- 
trix elements are real because of time-reversal symme- 
try. (This is the Gaussian orthogonal ensemble p^] .) A 
small opening in the cavity is described by a real, non- 
random M X N coupling matrix W, with N the number 
of wave channels transmitted through the opening. (For 
an opening of area A, N ~ 2nA/\ 2 at wavelength A.) 
Modes of the open cavity are complex eigenvalues (with 
negative imaginary part) of the non-Hermitian matrix 
H = Hq — iirWW T . The scattering matrix S at fre- 
quency uj is related to H by O] 



5=1- 2niW T {uj - H^W . 



(2) 



It is a unitary and symmetric, random N x N matrix, 
with poles at the eigenvalues of H . 

We now assume that the cavity is filled with a ho- 
mogeneous amplifying medium (amplification rate l/r a ). 
This adds a term i/2r a to the eigenvalues, shifting them 
upwards towards the real axis. The lasing mode is the 
eigenvalue fl — iT/2 closest to the real axis, and the laser 
threshold is reached when the decay rate T of this mode 
equals the amplification rate l/r a |14[ . Near the laser 
threshold we need to retain only the contribution from 
the lasing mode (say mode number I) to the scattering 
matrix (g), 
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where U is the matrix of eigenvectors of H . Because H 
is a real symmetric matrix, we can choose U such that 
U^ 1 = U T and write Eq. (||) in the form 



1 



S nm = <y n a m {uj - fl + iT/2 - i/2r a ) 1 , (4) 

where cr„ = (— 2Tti) 1 / 2 (W T U) n i is the complex coupling 
constant of the lasing mode I to the n-th wave channel. 
The Petermann factor K is given by 



1 N 



(5) 



The second equality follows from the definition of a n [[15| , 
and is the matrix analogon of Siegman's non-orthogonal 
mode expression 0]. The first equality follows from the 
definition of K as the factor multiplying the Schawlow- 
Townes linewidth JTs) . One verifies that K > 1 because 
(U^U)u > (U T U) U = 1. 

The relation (||) serves as the starting point for a 
calculation of the statistics of the Petermann factor in 
an ensemble of chaotic cavities. Here we restrict our- 
selves to the case iV = 1 of a single wave channel, leav- 
ing the multi-channel case for future investigation. For 
N = 1 the coupling matrix W reduces to a vector a = 
(Wh, W 2 i, . . • , Wmi)- Its magnitude \a\ 2 = (MA/tt 2 )w, 
where w 6 [0, 1] is related to the transmission probability 
T of the single wave channel by T — Aw (I + w)~ 2 . We 
assume a basis in which Hq is diagonal (eigenvalues u> q ). 

If the opening is much smaller than a wavelength, then 
a perturbation theory in a seems a natural starting point. 
To leading order one finds 



K = 1 + (27rai) 2 ^ 



qjtl 



{uji - U q ) 2 



(6) 



The frequency O and decay rate T of the lasing mode are 
given by and 2na 2 , respectively, to leading order in 
a. We seek the average (K)a,r of K for a given value 
of and r. The probability to find an eigenvalue at 
Lu q given that there is an eigenvalue at uji vanishes lin- 
early for small \ui q — as a consequence of eigenvalue 
repulsion constrained by time-reversal symmetry. Since 
the expression (J6j) for K diverges quadratically for small 
\uj q — u>i\, we conclude that (-K")n,r does not exist in per- 
turbation theory. This severely complicates the problem. 

We have succeeded in obtaining a finite answer for the 
average Petermann factor by starting from the exact re- 
lation 



UqlZl = L) q U q l - ilTCtq ^ Oi p Upl 



(7) 



between the complex eigenvalues z q of H and the real 
eigenvalues cu q of Hq. Distinguishing between q = I and 
q 7^ I, and defining d q — U q i/Uu, we obtain two recursion 
relations, 



Z{ = uji — iiraf — iirai a q d q , 

q^l 



Tra 



. j " ^q 
ld q = | Oil 

Zl - iOq 



(8b) 



The Petermann factor of the lasing mode / follows from 

-l 
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(9) 



We now use the fact that zi is the eigenvalue closest to 
the real axis. We may therefore assume that z; is close 
to the unperturbed value u>i and replace the denominator 
zi - u q m Eq. © by lui — Lo q . That decouples the two 
recursion relations, which may then be solved in closed 
form, 



id q = 



u>i — inaf (1 + it: A) 1 , 
na q ai . i 

- (1 + ITT A) 



We have defined A — J2 q ^i a g( w / — u q) 1 
rate of the lasing mode is 

T = -2Imz ; = 27ra 2 (l + n 2 A 2 )' 1 



(10a) 
(10b) 

The decay 
(11) 



Since the lasing mode is close to the real axis, we may 
linearise the expression (^|) for K with respect to T, 



K = 1 + 4^(Imrf 9 ) 2 = 1 



q7 tl 



{2ttT/A)B 

1 + TT 2 A 2 



(12) 



with B = A^a^wi - w q ) 2 . 

The conditional average of K at given T and Vt, can be 
written as the ratio of two unconditional averages, 

(K)n,r = 1 + (27rr/A)(S(l + t: 2 A 2 )- 1 Z) / (Z) , (13a) 
Z = 6{n-uJi)5(r-2ira 2 (l + Tr 2 A 2 )- 1 ) . (13b) 

In principle one should also require that the decay rates 
of modes q ^ / are bigger than T, but this extra con- 
dition becomes irrelevant for r — > 0. For M —* oo the 
distribution of a q is Gaussian oc exp(— ^a 2 Tr 2 /wA) [Q. 
The average of Z over ai yields a factor (1 + 7r 2 A 2 ) 1//2 , 



(K) a ,r = 1 + (27rr/A) 



(Bjl+^A 2 )- 1 / 2 ) 

((1 +7T 2 A 2 ) 1 /2) 



(14) 



where only the averages over a q and uj q (q ^ I) remain, 
at fixed u)i = Q. 

The problem is now reduced to a calculation of the 
joint probability distribution P(A,B). This is a techni- 
cal challenge, similar to the level curvature problem of 
random-matrix theory (T^Jl^]. The calculation will be 
presented elsewhere, here we only give the result: 



P(A,B) 



i nr^A 2 + w 2 



6 \ 2w B 7 / 2 



■ exp 



w 
2B 



ir 2 A 2 



w 



(15) 



Together with Eq. (|l4|) this gives the mean Petermann 
factor 
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we find the mean K, 
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FIG. 1. Average Petermann factor K for a chaotic cavity 
having an opening with transmission probability T. The av- 
erage is performed at fixed decay rate P of the lasing mode, 
assumed to be much smaller than the mean modal spacing 
A. The solid curve is the result ([L6|) in the presence of 
time-reversal symmetry, the dashed curve is the result ( |2c| ) 
for broken time-reversal symmetry. For small T, the solid 
curve diverges oc InT -1 while the dashed curve has the finite 
limit of it/3. For T — 1 both curves reach the value 2n/3. 
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(16) 



in terms of the ratio of two Meijer G-functions. We have 
plotted the result in Fig. [|, as a function of T = 4w(l + 
w)- 2 . 

The non-analytic dependence of the average K on T 
(and hence on the area of the opening |1| ) is a striking 
feature of our result. For T -C 1, the average reduces to 



(K)q r = H r- In — 

\ />M 6 A r 



(17) 



The non-analyticity results from the relatively weak 
eigenvalue repulsion in the presence of time-reversal 
symmetry. If time-reversal symmetry is broken by a 
magneto-optical effect (as in Refs. pC|^]] |), then the 
stronger quadratic repulsion is sufficient to overcome the 
lu~ 2 divergence of perturbation theory and the average 
K becomes an analytic function of T. For this case, we 
find instead of Eq. (O) the simpler expression 



(iOfi,r = l + (27rr/A)- 



(B) 



(1 + 7T 2 A 2 ) 

Using the joint probability distribution 



(18) 



, A U 2 A 2 +w 2 Y 
P(A,B)= K > exp 



w 
B 



ir 2 A 2 



(19) 
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r 

A3(l + w 2 ) ' 
It is equal to (K)a,r 



(20) 



shown dashed in Fig. [l| It is equal to (K)a,r = 1 
ivrTr/A for T < 1. 

So far we have concentrated on the average Petermann 
factor, but from Eqs. Q, @, and @ we can compute 
the entire probability distribution of K at fixed T. We 
define k = (K-1)A/TT. A simple result for P(k) follows 
forT=l, 



4^ 



k 1/2 exp(-7r/K) 



(21) 



7T 

1+ 2^ 



exp(-i7r/K) , kT<1. (22) 



P(«) 

and for T < 1, 
= 12^ 

As shown in Fig. ^, both distributions are very broad 
and asymmetric, with a long tail towards large k [^2| . 
The most probable (or modal) value of K — 1 ~ TT/A 
is parametrically smaller than the mean value ( |l7| ) for 
T< 1. 

To check our analytical results we have also done a 
numerical simulation of the random-matrix model, gen- 
erating a large number of random matrices Hq and com- 
puting K from Eq. (gb. As one can see from Fig. ||, the 
agreement with Eqs. ([2l|) and ( p2| ) is flawless. 

In conclusion, we have shown that chaotic scatter- 
ing causes large statistical fluctuations in the quantum- 
limited linewidth of a laser cavity. We have examined in 
detail the case that the coupling to the cavity is via a 
single wave channel, but our random-matrix model ap- 
plies more generally to coupling via an arbitrary number 

of wave channels. We have computed exactly the dis- 
tribution of the Petermann factor for N = 1. It remains 
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FIG. 2. Probability distribution of the rescaled Petermann 
factor k = (K - 1)A/PT for T = 1 and T < 1. The solid 
curves follow from Eqs. ( pi] ) and ( p2[ ) . The data points follow 
from a numerical simulation of the random-matrix model. 
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an open problem to do the same for N > 1. This prob- 
lem is related to several recent studies of the statistics 
of eigenfunctions of non-Hermitian Hamiltonians p3| , p4T] , 
but is complicated by the constraint that the correspond- 
ing eigenvalue is the closest to the real axis. Our study 
of a system with a fully chaotic phase space complements 
previous theoretical work on systems with an integrable 
dynamics. Chaotic laser cavities of recent experimental 
interest p5| have a phase space that includes both inte- 
grable and chaotic regions. The study of the quantum- 
limited linewidth of such mixed systems is a challenging 
problem for future research. 
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European Union. 



[1] A. L. Schawlow and C. H. Townes, Phys. Rev. 112, 1940 
(1958). 

[2] It is assumed that F is much less than the linewidth of 
the atomic transition and also that the lower level of the 
transition is unoccupied. 

[3] K. Petermann, IEEE J. Quantum Electron. 15, 566 
(1979). 

[4] A. E. Siegman, Phys. Rev. A 39, 1253 (1989); 

Phys. Rev. A 39, 1264 (1989). 
[5] Y.-J. Cheng, C. G. Fanning, and A. E. Siegman, 

Phys. Rev. Lett. 77, 627 (1996). 
[6] M. A. van Eijkelenborg, A. M. Lindberg, M. S. Thijssen, 

and J. P. Woerdman, Phys. Rev. Lett. 77, 4314 (1996). 
[7] M. Brunei, G. Ropars, A. Le Floch, and F. Bretenaker, 

Phys. Rev. A 55, 4563 (1997). 
[8] P. Grangier and J.-P. Poizat, Eur. Phys. J. D 1, 97 

(1998). 

[9] A. E. Siegman, Phys. Rev. A, to be published. 
[10] F. Haake, Quantum Signatures of Chaos (Springer, 
Berlin, 1991). 

[11] J. U. Nockel and A. D. Stone, Nature 385, 45 (1997). 
[12] M. Mehta, Random Matrices (Academic, New York, 
1990). 

[13] J. J. M. Verbaarschot, H. A. Weidenmiiller, and M. R. 
Zirnbauer, Phys. Rep. 129, 367 (1985). 

[14] T. S. Misirpashaev and C. W. J. Beenakker, Phys. Rev. A 
57, 2041 (1998). 

[15] To prove the second equality in Eq. (^), write ^2 n \o- n \ 2 = 
2Tr(U t WW T U)u = 2i{U\H-H )U)u = 2i{U^U)u{uo - 
iF/2) — 2i(W HoU)u and take the real part. 

[16] For the first equality in Eq. (|s|), write the linewidth 
Suj = T — l/r a in terms of the output power / = 
ftrSSUu/2n = (£ n \a n \ 2 ) 2 (T - l/ra)- 1 = KT 2 /5uj. 
The linewidth differs from the Schawlow- Townes value 
(|l|) by a factor 2K. The extra factor 2 arises from the 
suppression of amplitude fluctuations in the non-linear 



[17] 
[18] 
[19] 

[20] 

[21] 
[22] 

[23] 
[24] 
[25] 



regime above the laser threshold, as explained by P. Gold- 
berg, P. W. Milonni, and B. Sundaram, Phys. Rev. A 44, 
1969 (1991). 

F. von Oppen, Phys. Rev. Lett. 73, 798 (1994); 
Phys. Rev. E 51, 2647 (1995). 

Y. V. Fyodorov and H.-J. Sommers, Z. Phys. B 99, 123 
(1995). 

The transmission probability T is related to the area A of 
the opening by T ~ A 3 /X 6 for T < 1. See H. A. Bethe, 
Phys. Rev. 66, 163 (1944). 

H. Alt, H.-D. Graf, H. L. Harney, R. Hofferbert, H. 
Lengeler, A. Richter, P. Schardt, and H. A. Wei- 
denmiiller, Phys. Rev. Lett. 74, 62 (1995). 
U. Stoffregen, J. Stein, H.-J. Stockmann, M. Kus, and F. 
Haake, Phys. Rev. Lett. 74, 2666 (1995). 
For the case of broken time-reversal symmetry, we find 



P(k) = 



1^2 



-1/2-5/2 



exp(-27r//t) for T 

21 



, and P(k) = 

k "'"||+ 7t/2k + (tt/2k) j ~\ exp (— 7r/2«;) for 
T < 1. 

J. T. Chalker and B. Mehlig, Phys. Rev. Lett. 81, 3367 
(1998). 

R. A. Janik, W. No erenberg, M. A. Now ak, G. Papp, and 
I. Zahed, preprint ( |sond-mat/9902314| ). 
C. Gmachl, F. Capasso, E. E. Narimanov, J. U. Nockel, 
A. D. Stone, J. Faist, D. L. Sivco, and A. Y. Cho, Science 
280, 1556 (1998). 



4 



